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We study the parity effect and transport due to quasiparticles in circuits comprised of many su¬ 
perconducting islands. We develop a general approach and show that it is equivalent to previous 
methods for describing the parity effect in their more limited regimes of validity. As an example 
we study transport through linear arrays of Josephson junctions in the limit of negligible Josephson 
energy and observe the emergence of the parity effect with decreasing number of non-equilibrium 
quasiparticles. Due to the exponential increase in the number of relevant charge states with increas¬ 
ing length, in multi-junction arrays the parity effect manifests in qualitatively different ways to the 
two junction case. The role of charge disorder is also studied as this hides much of the parity physics 
which would otherwise be observed. Nonetheless, we see that the current through a multi-junction 
array at low bias is limited by the formation of meta-stable even-parity states. 

PACS numbers: 73.23.Hk,85.25.Cp,73.23.-b 


In superconducting circuits of small dimensions, charg¬ 
ing effects play an important role. On the one hand 
the Coulomb blockade leads to charge pinning and an 
effective suppression of electronic transport at low bias 
voltage. On the other hand the superconducting nature 
manifests in the parity effect, i.e. given a odd number of 
electrons on the superconductor there is one remaining 
quasiparticle dominating the transport properties in the 
low bias regimeip— . Strictly speaking this picture is true 
for equilibrium and very low temperatures. However if a 
non-equilibrium situation is imposed, e.g. by applying a 
finite bias voltage, the average number of quasiparticles 
may be increased. Recently such non-equilibrium quasi¬ 
particle effects have been investigated in superconducting 
qubits^i^— and single electron transistors (SETs)^“— . In 
this context the interplay of charge transport, the exci¬ 
tation of non-equilibrium quasiparticles and the observa¬ 
tion of the parity effect has been the subject of recent 
experiments with SETsi^. Based on related theoretical 
modellingiSr— we extend the prevailing transport theory 
of multi-junction circuits and show that this approach 
removes the ambiguities of previous approaches when in¬ 
cluding parity effects for more than one superconduct¬ 
ing island, although our approach is equivalent to earlier 
work in the appropriate limits. As an example, we per¬ 
form the first analysis of the parity effect in linear multi¬ 
junction arrays and make a number of predictions for 
the electronic transport signatures that can be identified 
with the parity effect in these systems. 


I. THE PARITY EFFECT IN 
MULTI-JUNCTION CIRCUITS. 


Conventional equilibrium quasiparticle theory states 
that transport through a Josephson junction is exponen¬ 



FIG. 1. (a) Quasiparticle tunnelling rate (Eq.[5} as a function 
of energy difference between charge configurations {5E) for 
four different parity configurations, ie. odd or even charge 
states on the source and destination islands respectively. The 
rates are evaluated at Te = 222 mK which is just below the 
parity temperature T* = 269 mK for these parameters - see 
text for details, (b) Parity rates evaluated in the middle of 
the sub-gap region [5E = A) as a function of Te. The scaling 
behaviour for the four rates as a function of NjkBTe. can be 
clearly seen above and below the parity temperature. 


tially suppressed when the applied bias V across the junc¬ 
tion is less than twice the superconducting gap, A. How¬ 
ever upon closer inspection there is a measurable sub¬ 
gap transport contribution for V < 2A, which depends 
strongly on the charge state parity of the islands either 
side of the junction. This parity effect is traditionally 
modelled in single-electron transistors using one of two 
approaches. 

In the first approach, the free energy of the circuit 
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gains an additional free energy contribution (which de¬ 
pends on temperature, superconducting gap and island 
volume) due to the parity of the occupancy of the is- 
Iand^i^iii^i22ri22. In the alternative approacbSii^i^i^^, the 
rates associated with both the equilibrium quasiparticles 
and the odd quasiparticle state ‘at gap’ must be com¬ 
puted separately taking into account the relative chem¬ 
ical potential differences between island and lead. Al¬ 
though in SETs both approaches describe similar physics; 
in the multi-junction case significant complications (both 
conceptual and technical) arise when applying either of 
these methods. In this work we show how to describe the 
parity contribution in a general way, as well as show how 
these earlier methods are considered limiting cases of the 
theory as presented here. 

Following more recent work on normal¬ 
superconducting-normal SETsi^, sub-gap quasiparticle 
effects can be included in a consistent way such that the 
rate for an arbitrary charge transfer event is computed 
based on the initial parity of the origin and destination 
islands. Throughout this discussion, we parameterize 
the distribution of non-equilibrium quasiparticles by an 
effective electron temperature Tg, rather than explicitly 
keeping track of the non-equilibrium quasiparticle num¬ 
ber on each island as was considered in Ref. [l^. This 
formalism is also applicable when modelling the parity 
corrections to sub-gap quasiparticle transport at base 
temperatures beyond the regime typically associated 
with the parity effect itself. 

To include the contribution from parity dependent 
quasiparticle tunnelling in a general way, we scale the 
Fermi function for each island in the circuit by a factor 
Aj, which depends on whether the charge state of the 
island n is odd or even (j = n mod 2). This accounts 
for the fact that in the odd charging state there is at 
least one quasiparticle remaining unpaired (see Ref. 
for further details). This modification of the Fermi func¬ 
tion f{E) —>■ Ajf{E), is strictly true only for E > 0 and 
we define. 


Aj = [tanh(Afgp)]^ (1) 
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This rate is expressed as a function of the energy dif¬ 
ference 6E between initial (n, m) and final charge states 
(n — 1,771 -I- 1 or 77 -I- 1, 777 — 1), where feiE) is defined 
as the Fermi function at temperature Tg and Rt is the 
junction normal tunnel resistance. In this context n and 
777 indicate the initial even-odd parity of the origin and 
destination islands respectively. 

When considering the movement of a single charge be¬ 
tween two islands, the rate given by Eq. [5] depends on 


where we parameterise the average number of excited 
quasiparticles Nqp a o^^d^ 


Nqp 
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A 
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Here, the average number of excited quasiparticles is ex¬ 
pressed in terms of the superconducting density of states 
evaluated at the Fermi level N{0), the volume of the 
island V, the superconducting gap A and the effective 
temperature Tg of the quasiparticles. When the num¬ 
ber of excited quasiparticles (Eq. is less than one, the 
scaling factor Aj shows markedly different behaviour for 
even and odd charging states, leading to parity depen¬ 
dent transport signatures. The crossover temperature 
T* below which these effects can be observed is given by, 

T* c- _^_ ('.Ti 

kB ln(iV(0)FV27rAfcB'r*) ^ ^ 

which must be solved self-consistently. In the limit of 
electron temperature > T*, the even-odd distinction 
vanishes and therefore the parity effect is unobservable. 

Although expressing the non-equilibrium contribution 
in terms of Aj is a very general approach, for f{E > 0) <C 
1 it proves to be both conceptually and computationally 
useful to parameterise the non-equilibrium quasiparticle 
distribution in terms of a modified chemical potential, 
fij. To do this we express Aj in the form of a shifted 
Fermi distribution such that 


Hj{E) = sgn{E){-iy kBTe In [tanh(7V,p)] (4) 

where the factor (—1)-^ takes into account the odd-even 
discrepancy and the sgn(i?) term accounts for the fact 
that the f{E) —>• Ajf{E) replacement applies strictly to 
positive energy differences. 

To study parity effects in a general way which will be 
applicable to multi-junction circuits, the single-electron 
tunnelling rate between islands also possess a even/odd 
charge state dependence and in general is given by 


fe{E - fin{E))[l - fe{E + SE - fim{E + ,5^))]. (5) 


the initial parity of both islands - giving four possible 
rates. If we consider the overall scaling of the sub-gap 
rates (insert to Fig.[T|), we see that above the parity tem¬ 
perature Tg > T*, all four rates scale oc exp[—A/fc^Tg]. 
Below the parity temperature, the ‘even-even’ rate scales 
oc exp[—2A/fesTg], whereas the other three rates are ap¬ 
proximately temperature independent for Tg < T*. 

To understand this behaviour, we can approximate 
Eq. [5] by expanding around the divergences in the BCS 














density of states, N{E)/N{0). Expanding E = A(1 + e) where 
and taking the dominant terms for each side of the den¬ 
sity of states, we obtain. 
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rn,m{^E) = 


A N{A + SE) 
e^Rr iV^O) 


[ de^g{E,SE) ( 6 ) 

Jo 


g{E, SE) — /(—A(1 -|- e) — SE -I- ^n)[l — + e) + 


We can then evaluate g{E, SE) and therefore the integral 
in various limits of interest. 

In the limit Tg > T*, the parity dependent chemical 
potential term « 0 and therefore g{E,SE) simplifies 
considerably. If we assume A ^ ksTe and SE > eA, we 
obtain g{E,SE) exp[—eA/ZcsTe] and therefore 

p(Te>T*) _ iV(A -I- SE) Ngp 

^ N{0) e’^RTN{Q)V ^ ’ 

The sub-gap quasiparticle rate as a function of SE there¬ 
fore takes on the functional form of the BCS density-of- 
states near the divergence, for temperatures above the 
parity temperature. The magnitude of this rate scales 
proportional to the quasiparticle number Ngp and as ex¬ 
pected is independent of the parity of the source and 
destination islands. 

Turning to the low temperature case (Tg T*), we 
now must deal with the various values of g.n,m- Taking 
the limit of Ngp <C 1 we can express the chemical poten¬ 
tial shift as 


(-l)"[fcBTgln(iVeff)-A] (9) 

where Wfi = N{Q)Vy/^UTKkBTl. Using this expression, 
we then evaluate g{E,SE) for each of the four cases, in 


the limit that A ^ fcsTe, SE > eA and Ngs ^ 1. 
gives the following expressions. 

This 

p(re<T*) ^ 

N{A + SE) Nip 

iV(0) e'^RTN{Q)V 

(10) 

p(Te<T*) ^ 
^ 0,0 

_N{A + SE) 1 

iV(0) e^RTN{0)V 

(11) 

p(Te<T*) p(Te<T*) 

p(Te<T*) p(Te<T*) ^ 0,0 -l-le.e 

^ o,e ^ e,o 2 

(12) 

In all rates, we see 

the characteristic density of states 


dependence on SE as well as a factor of N^p difference 
between the odd-odd and even-even rates. Therefore 
the even-even rate scales with Ngp oc exp(—2A/fcBTe) 
whereas the odd-odd rate is approximately constant ap¬ 
proaching zero temperature. As can be seen in Fig.[TJ the 

(T <^T *) 

odd-even and even-odd rates are equal and ~ ro,o /2 


Mm)] + /(A(l + e) — Mn)[l ~ /(^(l + c) + SE — g,m)]- (7) 


(as the contribution from the even-even rate is negligible 
at low temperatures). 

The approximate rates given above can be compared to 
previous work on the parity effect in simpler circuits. The 
mapping to the free energy shift of A often ascribed to 
the odd charge stat o^di^^~ 2i follows directly from Eq. O 
However, we can immediately see the lack of generality 
of that approach because the sign of the shift depends 
on the parity of both the source and destination charge 
states. In our case of many islands a simple A correction 
for each odd charge state is manifestly not sufficient. 

Previous work on applying a shifted chemical poten- 
tial2S“— is in principle similar to our approach. In that 
case the “odd to even” transition rate (meaning the tran¬ 
sition from an odd to even state of the same island) is 
equivalent to our Eq. 1111 However it is not clear how 
to easily generalise this method to multiple islands when 
the parity of both the source and destination islands must 
be taken into account. Furthermore, Eos. ITOlfT^ are only 
approximations to the general expression Eq. O due to 
the relatively crude approximation to the integral over 
e. This becomes particularly important when comparing 
quantitatively to experiment. 


II. THE PARITY EFFECT IN JOSEPHSON 
JUNCTION ARRAYS 

A Josephson junction array (JJA) is the multi-junction 
generalisation of the (superconducting) single-electron 
transistor. Increasing the number of junctions changes 
the electrical response of such a circuit markedly^ir— 
when compared to the simple two junction devices. New 
and interesting effects are observed, including hystere¬ 
sis^, soliton propagation2ii2^, non-trivial magnetic field 
effects^i^i^ and correlated electron transpor t^^i^^'^^ . 
Detailed understanding of junction arrays also promise 
new superconducting devices; such as qubit designs 
based on large kinetic inductance^^— , or terahertz ra¬ 
diation sources^^—. Although the qualitative theory of 
Josephson junction arrays has been established for some 
time^i^ij direct quantitative comparison between the¬ 
ory and experiment in these devices is still elusive and 
fraught with difficulty as they display qualitatively dif¬ 
ferent physics than that seen in single-electron transistors 
(SETs) and other few junction devices. 
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(a) Cj,Rt 



FIG. 2. (a) JJA circuit under consideration, consisting of a 
linear chain of Josephson junctions with Josephson energy Ej 
and capacitance Cj. The circuit is driven by a symmetrically 
applied voltage source V and each junction sees an effective 
capacitance to ground Cg- (b) Current through the array 
(logarithmic scale) as a function of applied voltage. Effective 
electron temperature is varied from below to above the parity 
crossover temperature, which is T* = 269 mK for these pa¬ 
rameters. The lines in the figure are the average of 10 KMC 
runs, each of which consists of 10® events. The maximum and 
minimum of these 10 runs are indicated by the pale shading. 


To illustrate the manifestation of the parity effect in 
multi-island circuits, we simulate transport through a lin¬ 
ear JJA, see Figl^Ka). We employ the kinetic Monte- 
Carlo metho d^^'^^~ — following the procedure as detailed 
in Refill. The energy of various charge states of the array 
is computed based on purely electrostatic considerations 
as we confine our investigation to the negligible Ej limit 
and therefore do not consider Cooper-pair transport (as 
Rt = 1M12 ^ Rq = hje'^). We assume a supercon¬ 
ducting gap of A = 200/xeV, normal density of states at 
the Fermi energyi^*^ -^(0) = 1.4477 x 10“^^ m“^ J“^ and 
island volume V = 0.0014/im^, consistent with experi¬ 
ments on aluminium based JJA a^^i^^ . Throughout this 
discussion, we consider a JJA of length iV = 50 with 
a ratio of junction capacitance Cj = 0.5fF to ground 
capaci tance C g = 20aF that gives a soliton length^i^i 
A = y/Cjj Cg = 5. In this regime the array can be con¬ 
sidered ‘long’ although correlated transport effects are 
still importan t^°i^^ . 

The ‘smoking gun’ of parity effects in superconduct¬ 
ing SETs is the observation of current plateaus at low 
temperature for odd charge states, while the current is 
suppressed completely for even charge states^ii^. We be¬ 
gin by considering the equivalent experiment for a JJA. 
Fig. [5] shows the I-V characteristics for a JJA as a func- 



FIG. 3. (a) Current as a function of voltage and electron tem¬ 
perature for maximal disorder. Below the parity temperature, 
the current is strongly suppressed at low voltages. The oscil¬ 
lations in the current due to parity effects are also lost due 
to disorder, (b) The current at Te = 170 mK as a function 
of disorder strength. In this subfigure we ensemble averaging 
over 50 different disorder realisations, the variance of which is 
show via pale shading. Although there is considerable varia¬ 
tion with disorder realisation, qualitatively we see that using 
t; = 1 insures numerical convergence for a given realisation. 


tion of effective electron temperature. For Te < 170 mK, 
we see only sporadic (or no) conduction as the system 
is too easily trapped in meta-stable states. As Te in¬ 
creases, the characteristic suppression of current due to 
the parity effect manifests as oscillations in the current 
at low bias. These oscillations stem from the interplay 
between the parity dependent tunnelling rates and the 
voltage dependent filling factors on each island. When 
V > 22mV we see the step in current associated with 
the breaking of Cooper-pairs at every junction, ie. when 
V>Nx2A = 20mV. 

Interestingly, the parity oscillations vanish for V > 
13mV although the magnitude of the current is still ap¬ 
proximately constant below T* but increases rapidly for 
Te > T*. At hrst glance this would appear to be a tran¬ 
sition associated with V = N x A however it actually 
depends on the interplay of charging energy and parity 
effects. At low voltages, dipole states can form which are 
stable for certain combinations of voltage parity and it 
is these states which block the flow of current. Above 
a certain voltage, these metastable states can dissociate 
via interactions with neighbouring charges - leading to 
more robust conduction. One can think of this in terms 
of a phase-space argument where the number of available 
states grows with increasing voltage, therefore allowing 
the system to avoid getting trapped in local minima. 

Experimentally, background charges within the device 
and substrate lead to random offset charges^^— . We 
can model this disorder as random offset charges on the 
islands of the JJAi^, \qbk\/e < rj, which we assume to 
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FIG. 4. (a) Current plotted as a function of Te for various val¬ 
ues of V. Here we see at low voltages a characteristic crossover 
from 2A/kBTe (dashed-dotted line) to A/ksT^ (dashed line) 
scaling, indicating parity limited conduction. At higher volt¬ 
ages, the increased state space means that the current is no 
longer limited by the even-even rate but is instead dominated 
by the odd parity processes, showing an initial constant scal¬ 
ing with temperature. Irrespective of voltage (while still in 
the sub-gap region), at temperatures above T* the current 
scales with the parity independent rate, (b) The inset shows 
the same data plotted as a function of inverse temperature, 
which clearly shows the two regimes. 


be static on the time scale required to measure a single 
current point. In Fig. [3] we see the effects of increasing 
disorder is to suppress conduction for small bias, as well 
as eliminating the parity dependent oscillations as a func¬ 
tion of voltage bias. We see convergence of the response 
as a function of disorder strength (inset to Fig. [S]), with 
V ^ 0 -'^ being sufficient to model ‘maximal disorder—. 
As well as the issues of rare-events at low bias discussed 
for the non-disordered case, in the inset to Fig. [3] we com¬ 
pute the variance in the current as a function of 50 dif¬ 
ferent disorder realisations. We see the same qualitative 
suppression of the parity oscillations with increasing dis¬ 
order, even taking into account the variance due to disor¬ 
der. For all subsequent simulations (as well as the main 
Fig. [31) we take the disorder realisation to be constant 
(with 77 = 1) for a given I-V curve. This corresponds to 
the case where the charge disorder is maximal but stable 
over the entire experimental timescale of interest. 

Although the characteristic parity oscillations are not 
visible in the presence of maximal disorder, the overall 
scaling behaviour of the current is still a strong function 
of the parity dependent rates. This scaling behaviour 
at low bias can be simply understood in terms of the 
scaling of the NEQ rates shown in Fig. [TJ To observe 
the scaling of the tunnelling rate above and below T* ie. 
oc exp(—A/fegTe) and oc exp{—2A/kBTe) respectively; 
in FiglHwe plot the current at fixed voltage. In the low 
bias regime {V = 5mV) we see clear evidence of the cross¬ 
over from 2A/fcTe to A/kT^ scaling. This regime is where 



FIG. 5. (a) Fourier transform of the charge-charge correlation 
function as a function of temperature. As the temperature is 
increased, the mean transition rate and therefore the current 
also increases. As the temperature is increased, correlated 
conduction sets in soon after conduction itself begins. The 
ratio of correlation frequency to current /peak 11 = e is con¬ 
sistent with the charge carriers being single electrons, (b) 
The width of the correlation peaks normalised by the posi¬ 
tion (the reduced line width) shows a particularly clear signa¬ 
ture of the parity effect. As a function of temperature we see 
the reduced line width plateau above the parity temperature 
(indicated with triangles for different values of the supercon¬ 
ducting gap). 


charge-charge correlation effects are strongest due to the 
low filling of the array^^i^i^ which is also the regime 
that should show parity effects most clearly. 

In the low bias regime, we interpret this cross-over as 
single isolated charges moving through the array largely 
independently. However, due to the background charge, 
the effective parity of the islands they encounter varies 
and therefore the slowest rate (the even-even rate) is the 
limiting factor. The scaling behaviour therefore mimics 
that of the even-even rate below and above the parity 
temperature. At higher bias, the transition from con¬ 
stant to A/kTf, scaling illustrates that the dominant tun¬ 
nelling events are associated with odd-odd and odd-even 
processes. In this case the higher applied bias pushes 
charges closer together which increases the available state 
space and the system can escape from even parity states 
which would otherwise limit conduction. Such a transi¬ 
tion from constant to A/kTe scaling was recently seen^ 
in zero-bias conductance experiments on flux tuneable 
Josephson junction arrays - which suggests a thermal 
quasiparticle origin for the temperature dependence of 
the conductance. 

Recent experiments have demonstrated counting of in¬ 
dividual electrons within a - which raises the 

question of whether signatures of parity can be seen in 
the charge-charge correlations during transport. Even in 
conventional (normal) conducting JJAs, there is strong 
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interplay between applied voltage and correlated trans¬ 
port (through the average charge density). To focus on 
the role of parity, we compute the charge-charge corre¬ 
lation functio n^^i^° on site number 15 of the array for a 
fixed applied voltage V = 2mV, sampled with a band¬ 
width of 5 MHz. As a function of temperature (at fixed 
voltage) we see three distinct regions, see Fig. 

At very low temperatures and currents, we see no cor¬ 
relation as the transport is too slow on the time scale of 
the simulations. As the temperature and therefore cur¬ 
rent increases, strongly correlated transport sets in with 
the correlation peak frequency scaling linearly with cur¬ 
rent according to I = e/peak, reflecting the fact that the 
charge carriers are single electrons. Ultimately the am¬ 
plitude of the correlation peak reduces due to increasing 
charge noise at high currents - showing a surprising simi¬ 
larity as a function of temperature to the ‘washout’ seen 
at high voltage bias^. 

More subtly, the role of the parity effect can be seen 
in the ‘reduced line width’ of the response, ie. peak 
width/peak position, as a function of temperature (in¬ 
set to Figl5|). Once correlated transport sets in, the peak 
width reduces as a function of increasing temperature 
until it reaches a constant value when Ri T*, which is 
largely independent of the value of the superconducting 
gap. We can ascribe this step behaviour to the addi¬ 
tional noise in the correlation signal due to the parity 
effect, which effectively vanishes once Tg reaches T*. 


III. CONCLUSION 

Josephson junction arrays provide a tantalising play¬ 
ground for studying many-body effects as they are a con¬ 
trollable, artificial system which is truly one-dimension 
and yet displays correlation electron effects. It is there¬ 
fore supremely disappointing that experimental results 
to date can only be explained qualitatively at best. An 
important contributor to this situation is the difficulty in 
(experimentally) filtering out or (theoretically) account¬ 
ing for quasiparticle effects. We have shown how to model 
the transition from many to single quasiparticle excita¬ 
tions in multi-junction circuits, in particular we consid¬ 
ered the quasiparticle parity in a consistent way. Even 
when considering the limit of strongly disordered offset 
charges, the crossover to the parity regime can be ob¬ 
served in both the current-voltage characteristics and the 
charge-charge correlation function. 
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